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Abstract 

The 0{a) radiative corrections to the process 7*, Z (l = e, fi) 

are calculated. Factorizing the collinear singularity associated with initial 
state photon bremsstrahlung into the parton distribution functions, we find 
that initial state corrections have a much smaller efi'ect than final state ra- 
diative corrections. Due to mass singular logarithmic terms associated with 
photons emitted collinear with one of the final state leptons, QED radiative 
corrections strongly affect the shape of the di-lepton invariant mass distribu- 
tion, the lepton transverse momentum spectrum, and the forward backward 
asymmetry, Aps- They lead to a sizeable shift in the Z boson mass extracted 
from data, decrease the di-lepton cross section by up to 10%, and increase the 
integrated forward backward asymmetry in the Z peak region by about 7% at 
the Tevatron. We also investigate how experimental lepton identification re- 
quirements modify the effect of the QED corrections, and study the prospects 
for a high precision measurement of sin^ ^^^^ using the forward backward 
asymmetry at the Large Hadron Collider (LHC). 
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I. INTRODUCTION 



Over the last few years, the Standard Model (SM) of electroweak interactions has been 
successfully tested at the one-loop level. Experiments at LEP and the SLC have de- 
termined the properties of the Z boson with a precision of 0.1% or better, and correctly 
predicted the range of the top quark mass from loop corrections [0]. Currently, the Z boson 
mass is known to ±2.0 MeV, whereas the uncertainty of the W mass, Mvi/, is ±80 MeV 0. 
A precise measurement of Mv^ and the top quark mass, mtop, would make it possible to 
derive indirect constraints on the Higgs boson mass, M/^, via top quark and Higgs boson 
electroweak radiative corrections to Mvi/ [^. With a precision of 30 MeV (10 MeV) for the 
W mass, and 2 GeV for the top quark mass, Mu can be predicted with an uncertainty of 
about 50% (20%) IQ. Comparison of these constraints on Mn with the mass obtained from 
direct observation of the Higgs boson in future collider experiments will be an important 
test of the SM. 

A significant improvement in the W mass uncertainty is expected in the near future from 
measurements at LEP II and the Fermilab Tevatron collider The ultimate precision 
expected for from the combined LEP II experiments is approximately 40 MeV 0|. At 
the Tevatron, integrated luminosities of order 1 fb~^ are envisioned in the Main Injector 
Era, and one expects to measure the W mass with a precision of approximately 50 MeV |^ 
per experiment. The prospects for a precise measurement of M^y would further improve if 
a significant upgrade in luminosity beyond the goal of the Main Injector could be realized. 
With recent advances in accelerator technology Tevatron collider luminosities of order 
10^'^ cm~^ s~^ may become a reality, resulting in integrated luminosities of up to 10 fb^^ per 
year. With a total integrated luminosity of 30 fb~^, one can target a precision of the W 
mass of 15 - 20 MeV 0]. A similar or better accuracy may also be reached at the Large 
Hadron Colhder (LHC) 0. 

The determination of the W mass in a hadron collider environment requires a simulta- 
neous precision measurement of the Z boson mass, M^, and width, V z- When compared to 



2 



the value measured at LEP, the two quantities help to accurately determine the energy scale 
and resolution of the electromagnetic calorimeter, and to constrain the muon momentum 
resolution PQ]. 

Analogous to the W mass, a very high precision measurement of the effective weak mixing 
angle, sm^ O^^Jj |]TD[, can be used to extract information on the Higgs boson mass P,P 



At hadron colliders, the effective weak mixing angle can be determined from the forward 
backward asymmetry, ApB-, in di-lepton production in the vicinity of the Z pole [[T2| . 

In order to measure ApB and the Z boson mass with high precision in a hadron col- 
lider environment, it is necessary to fully understand and control higher order QCD and 
electroweak corrections. A complete calculation of the full 0{a) radiative corrections to 
P^P ^7*) Z ^ has not been carried out yet. In a previous calculation, only the final 



state photonic corrections had been included using an approximation in which the 



sum of the soft and virtual part is indirectly estimated from the inclusive Z — > £+£"(7) 

width and the hard photon bremsstrahlung contribution. 

In this paper, we present a more complete calculation of the 0{a) QED corrections to 
pj) 'J*, Z ^ £^£^. Real and virtual initial and final state corrections, as well as the 
interference between initial and final state corrections are included. Purely weak corrections 
are expected to be very small and are therefore ignored. Our calculation also takes into 
account the mass of the final state leptons, which regularizes the coUinear singularity associ- 
ated with final state photon radiation. Both Z and photon exchange diagrams with all 7 — Z 
interference effects are incorporated. The di-lepton invariant mass thus is not restricted to 
the Z peak region. Low mass Drell-Yan production is of interest because of the sensitivity 



to parton distribution functions (PDF's) at small x values |T^. High mass lepton pairs and 



the forward backward asymmetry above the Z peak can be used to search for additional 



neutral vector bosons, and to constrain their couplings |T^,|T8|. Results from our calculation 
have been used in Ref. [|T8| to compare experimental data with the SM prediction for ApB- 
To perform our calculation, we use the Monte Carlo method for next-to-leading-order 



(NLO) calculations described in Ref. |jT9|. The matrix elements for radiative Z production 



and decay are taken from Ref. and With the Monte Carlo method, it is easy to 
calculate a variety of observables simultaneously and to simulate detector response. Special 
care has to be taken in calculating the radiative corrections associated with photon radiation 
from the incoming quarks and antiquarks. In the parton model, quarks are assumed to be 
massless, and initial state photon radiation results in collinear singularities. The singular 
terms are universal to all orders in perturbation theory and can be removed by universal 



collinear counterterms generated by 'renormalizing' the parton distribution functions |l22| , |23 
in complete analogy to gluon emission in QCD. A calculation of QED corrections using def- 
inite, non-zero, values for quark masses and not factorizing the corresponding collinear 



logarithms leads to a considerable overestimation of the effects of initial state photon cor- 
rections. However, QED corrections to the evolution of the parton distribution functions 
are not included in our calculation; a complete fit of the PDF's including all QED effects 
is beyond the scope of this paper. The technical details of our calculation are described in 
Sec. II. 

Numerical results for pp collisions at ^/s = 1.8 TeV are presented in Sec. III. Due to 
the mass singular logarithms associated with final state photon bremsstrahlung in the limit 
where the photon is emitted collinear with one of the charged leptons, the di-lepton invariant 
mass distribution is strongly affected by QED corrections, in particular in the vicinity of 
the Z boson resonance. As a result, the value extracted for Mz from data is shifted to 
a lower value. The amount of the shift depends on the lepton mass, and the detector 
resolution QED radiative corrections also significantly affect the Z boson production 

cross section when cuts are imposed, the transverse momentum distribution of the leptons, 
and the forward backward asymmetry below the Z pole. For di-lepton masses between 
50 GeV and 100 GeV, the final state 0{a) QED corrections are larger than the 0{as) QCD 
corrections. 

In Sec. Ill, using a simplified model of the CDF detector as an example, we also inves- 
tigate how the finite energy and momentum resolution of realistic detectors affect the QED 
corrections. Electrons and photons which are almost collinear are difficult to discriminate, 



and the momenta of the two particles are thus recombined into an effective electron momen- 
tum 0,^] if they traverse the same calorimeter cell, or, alternatively, if their separation in 
the pseudorapidity - azimuthal angle plane is below a critical value. The second procedure 
completely eliminates the mass singular logarithms. With the first method, residual effects 
of these terms remain when both particles are almost coUinear, but hit different calorimeter 
cells. In practice, the numerical difference between the two procedures is moderate; in both 
cases the significance of the QED corrections is considerably reduced. In contrast, photons 
which are almost collinear with muons are rejected if they are too energetic |^ which results 
in residual logarithmic corrections to observable quantities in fi^fi^ production. Transverse 
momentum and rapidity cuts are found to affect the lepton pair invariant mass distribution 
and forward backward asymmetry in a similar way at the Born level and at 0{a^). 

Recently, it has been suggested that an ultra precise measurement of sin^ 6*1,^^ may 



be possible at the LHC {pp collisions at y/s = 14 TeV [^) in the muon channel, using the 
forward backward asymmetry in the Z peak region. At the LHC, the forward backward 
asymmetry is significantly reduced compared to the Tevatron because of the larger sea - sea 
quark parton flux. We find that the sensitivity of ApB to the effective weak mixing angle 
strongly depends on the rapidity range over which the leptons can be detected. The forward 
backward asymmetry at the LHC, including 0{a) QED and 0{as) QCD corrections, is 
studied in detail in Sec. IV. Finally, our conclusions are presented in Sec. V. 

II. METHOD OF CALCULATION 

The calculation presented here employs a combination of analytic and Monte Carlo 
integration techniques. Details of the method can be found in Ref. [|19|. The calculation of 



di-lepton production in hadronic collisions at 0{a^) includes contributions from the square 
of the Born graphs, the interference between the Born diagrams and the virtual one loop 



graphs, and the square of the real emission diagrams which we adopt from Refs. |20,21|. The 



diagrams contributing to the 0{a) QED corrections can be separated into gauge invariant 



5 



subsets corresponding to initial and final state corrections. The squared matrix element for 
the real emission diagrams is then given by 

|A<2^3|2 ^ |_^2^3|2 ^ 2Re[M'rHMy'r] + \My'\'. (1) 

A4]^^ and Ai^^^ are the separately gauge invariant matrix elements associated with initial 
and final state radiation. 

The basic idea of the method employed here is to isolate the soft and collinear singular- 
ities associated with the real photon emission subprocesses by partitioning phase space into 
soft, collinear, and finite regions. This is done by introducing theoretical soft and collinear 
cutoff parameters, 6s and 6c- Using dimensional regularization p6|, the soft and collinear 



singularities are exposed as poles in e (the number of space-time dimensions is = 4 — 2e 
with e a small number). In the soft and collinear regions the cross section is proportional 
to the Born cross section. The soft region is defined by requiring that the photon energy in 
the qq center of mass frame, E^, is < 6sV^/2 {s denotes the squared parton center of 
mass energy). We can then evaluate, in N dimensions, the 2^3 diagrams using the soft 
photon approximation, where the photon momentum is set to zero in the numerator, and 
integrate over the soft region. The soft singularities originating from final state photon radi- 
ation cancel against the corresponding singularities from the interference of Born and final 
state virtual corrections. Similarly, the soft singularities associated with initial state photon 
emission and interference effects between initial and final state radiation cancel against the 
corresponding singularities originating from initial state vertex corrections, and the Z'j and 
77 box diagrams, respectively. The remainder is then evaluated via Monte Carlo integration 
as part of the 2^2 contribution. For E^ > 6s\/S/2, the real photon emission diagrams are 
calculated in four dimensions p0|j21| using standard three body phase space Monte Carlo 



integration techniques. 

The collinear singularity associated with photon radiation from the final state lepton 
line is regulated by the finite lepton mass. The collinear singularities originating from 
initial state photon bremsstrahlung are universal to all orders of perturbation theory and 
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can be cancelled by universal collinear counterterms generated by renormalizing the parton 
distribution functions [^,^], in complete analogy to gluon emission in QCD ||2^. They 
occur when the final state photon and the partons in the initial state are collinear so that 
denominators of propagators such as 

i=-2p^.p^ (2) 

and 

u = -2pq ■ (3) 

vanish. Here, Pg {pq) denotes the quark (anti-quark), and p^ the photon four momentum 
vector. Only is divergent in the collinear limit; the initial - final state interference 

term, Re[M.'j~*^ {Miy^^)*] exhibits only soft singularities for massive final state leptons. In 
the collinear region, 1^1 < 5cS, |A^?^^p is evaluated in the leading pole approximation. 
After A^-dimensional integration over the photon phase space variables, the explicit singu- 
larity can be factorized into the parton distribution functions. The remainder is evaluated 
as part of the 2 — 2 contribution. If \u\ > ScS, the 2 — 3 diagrams are again evaluated 
numerically in four dimensions using the full three body phase space. 

In order to treat the 0{a) initial state QED corrections to di-lepton production in 
hadronic collisions in a consistent way, QED corrections should be incorporated in the global 
fitting of the PDF's. Current fits to the PDF's do not include QED corrections. A 



study of the effect of QED corrections on the evolution of the parton distribution functions 
indicates |2^ that the modification of the PDF's is small. We have not attempted to include 



QED corrections to the PDF evolution in the calculation presented here. The missing QED 
corrections to the PDF introduce an uncertainty which, however, probably is much smaller 
than the present uncertainties on the parton distribution functions. 

Absorbing the collinear singularity into the PDF's introduces a QED factorization scheme 
dependence. The squared matrix elements for different QED factorization schemes differ by 
the finite 0{a) terms which are absorbed into the PDF's in addition to the singular terms. 



As long as QED corrections to the PDF evolution are not included, the 0{a^) cross section 
will depend on the QED factorization scheme used. We have performed our calculation in 
the QED MS and DIS schemes, which are defined analogously to the usual MS ||2^ and 
DIS schemes used in QCD calculations. Unless noted otherwise, we will use the QED 
DIS scheme which minimizes the effect of the 0{a) QED corrections on the PDF by requiring 
the same expression for the leading and next-to-leading order structure function F2 in deep 
inelastic scattering. 

The 2 — > 2 contribution associated with initial state radiative (ISR) corrections, including 
the correction terms originating from the absorption of the initial state collinear singularity, 
can be obtained from the corresponding 0{as) QCD corrections by replacing {A/3)as 
by aQ"^, where Qg is the electric charge of the quark in units of the proton charge, in all 
relevant matrix element and cross section formulae. The 2 — 2 contribution induced by the 
soft and virtual final state radiative (FSR) corrections is given by: 
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where is the lepton mass and 
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Born 
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is the Born qq ^ 'y*, Z ^ i^i matrix element. Finally, the 2-^2 contribution induced by 
the 0{a^) initial - final state interference correction terms is given by 
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+ finite 77 and 7Z box terms 



(6) 



with 



ft,.. = log(|L) (7) 

and 

iz = — . (8) 

^ Mz ^ ^ 

In our calculation, we use the full s dependent width in the Z boson propagator. The ti 
and Ui are Mandelstam variables of the 2 — 2 reaction: 

ii = -'2pq-pe+, (9) 

ui = -2pg-pi-. (10) 

The finite terms from the 77 and jZ box diagrams are identical to those in e~^e~ qq and 



can be found in Refs. ||33| and |p3 



The end result of the calculation consists of two sets of weighted events corresponding 
to the 2^2 and 2-^3 contributions. Each set depends on the parameters Ss and 6c- The 
sum of the two contributions, however, must be independent of 6s and 6c, as long as the two 
parameters are taken small enough so that the approximations used are valid. In Figs. |l| 
and H we show the dependence of the pp — > £+£"(7) cross section in the Z peak region 
(75 GeV < m(£+£^) < 105 GeV) on 6s and 6c; 'm{i^i^) denotes the di-lepton invariant 
mass. To compute the cross section, we use here and in all subsequent figures the MRSA 



set of parton distribution functions , and take the renormalization scale fx and the QED 
and QCD factorization scales, Mqed and Mqcd, to be /i^ = Mq^j^ = Mq^jj-, = s. 

Figure |1| displays the cross section as a function of 6c (Fig. ^) and 6s (Fig. ^) for initial 
state radiative corrections only. In order to exhibit the independence of the cross section 
from the parameters 6s and 6c more clearly, we have not included the Born cross section in 
the 2 — s> 2 contribution. The ISR corrections to the cross section for electron and muon final 
states are virtually identical. While the separate 2^2 and 2^3 0{a) contributions vary 
strongly with 6s and 6c, the sum is independent of the two parameters within the accuracy of 
the Monte Carlo integration. The total contribution of initial state radiation diagrams to the 



total cross section in the Z pole region is found to be about 0.43% of the Born cross section 
for the parameters chosen. In the QED MS scheme, the contribution of the ISR diagrams 
is about 10% smaller than in the QED DIS scheme. QED corrections to the PDF's and 
purely weak one loop corrections to the matrix elements, both which are not included in our 
calculation, are expected to be of the same order of magnitude. 

In Fig. 1^, we show the •p'p £+£"(7) cross section in the Z peak region (75 GeV < 
m < 105 GeV) as a function of the soft cutoff parameter bs for electron and muon 
final states for FSR corrections. Radiation of photons coUinear with one of the leptons gives 
rise to terms proportional to \o<g[slvn\) log(5s) (see Eq. (^)) in both the 2^2 and 2 — > 3 
contributions. As demonstrated in Fig. ^, these terms cancel and the total cross section 
is independent of 5^. Due to the smaller mass of the electron, the variation of the 2—^2 
and 2^3 contributions with 5^ is more pronounced in the electron case. The solid line in 
Fig. indicates the cross section in the Born approximation. The total Oip?^ cross section 
in the e"'"e~(7) {fi^^^ {'-/)) case is found to be about 7% (3%) smaller than the Born cross 
section. The difference in the NLO 6+6^(7) and /i^yU^(7) cross section can be traced to 
residual logarithmic correction terms which arise from the finite lepton pair invariant mass 
range considered in Fig. |^ (see Sec. IIIA). If the integration would be carried out over the 
full range m(£+£^) > 2m£, these terms would vanish ^5|. From Fig. |^ one also observes 



that, due to the residual logarithmic terms, final state radiation effects are much larger than 
those which originate from initial state radiation. The 2—^2 and the 2^3 contributions 
to the FSR corrections each are trivially independent of the collinear cutoff 6c- 

Similar to the FSR corrections, one can show that the sum of the 2 — 2 and 2 
3 contributions of the initial - final state interference terms is independent of 6s- The 
interference terms are typically of the same size as the initial state corrections. 

The missing QED corrections to the PDF's create a dependence of the 0{a) initial state 
corrections on the factorization scale Mqed which is stronger than that of lowest order 
calculation. On the other hand, final state and initial - final state interference terms depend 
on the factorization scale only through the PDF's. These terms therefore exhibit a sensitivity 
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to the factorization scale which is similar to that of the lowest order calculation. Since the 
Born cross section and final state corrections are much larger than corrections from initial 
state radiation, the scale dependence of the complete 0{a^) cross section is similar to that 
of the Born cross section. 

In conclusion, in Figs. |I| and ^we demonstrated that the pp — > 7*, Z ^ £+£^(7), cross 
section for 75 GeV < m(£+£^) < 105 GeV is independent of the soft and collinear cutoff 
parameters 6s and 6c within the accuracy of the Monte Carlo integration. Independence 
of the cross section from these two parameters can also be demonstrated for lepton pair 
invariant masses below (m(£+£~) < 75 GeV) and above (m(£+£^) > 105 GeV) the Z peak. 
In the following, the soft and collinear cutoff parameters will be fixed to 6s = 10^^ and 
6c = 10~^, unless explicitly stated otherwise. 

III. 0{a) CORRECTIONS TO DI-LEPTON PRODUCTION AT THE TEVATRON 

We shall now discuss the phenomenological implications of 0{a) QED corrections to 
di-lepton production at the Tevatron {pp collisions at y/s = 1.8 TeV). We first discuss the 
impact of QED corrections on the lepton pair invariant mass distribution and the forward 
backward asymmetry. We then consider how the finite resolution of detectors and exper- 
imental lepton identification requirements modify the effects of the QED corrections, and 
investigate how 0{a) QED corrections affect the measured di-lepton {Z boson) cross section 
within the cuts imposed. Finally, we study the effect of the full radiative corrections on the 
Z boson mass extracted from data. The SM parameters used in our numerical simulations 
are Mz = 91.187 GeV, a(M|) = 1/128, Tz = 2.50 GeV and sin^ ^' J^* = 0.2319. These 
values are consistent with recent measurements at LEP, SLC and the Tevatron [Q]. 

A. QED Corrections to the Di-lepton Invariant Mass Distribution and ApB 

As we pointed out in Sec. II, final state photon radiation leads to corrections which 
are proportional to a\og{s/mj). These terms are large, and are expected to significantly 
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influence the shape of the di-lepton invariant mass distribution. The 0{a^) invariant 
mass distribution in the vicinity of the Z peak for the electron (solid line) and muon case 
(dotted line) is shown in Fig. ^together with the lowest order prediction (dashed line). No 
detector resolution effects or acceptance cuts are taken into account in any of the figures 
shown in this subsection. QED corrections decrease (increase) the cross section at (below) 
the peak. At the peak position, the differential cross section is reduced by a factor |^ 



with 

-K \ mf ) 

i.e. by about 30% in the electron case and by about 20% in the muon case. The shape of the 
Z boson resonance curve is seen to be considerably distorted by the 0{a) QED corrections. 
Photon radiation from one of the leptons lowers the di-lepton invariant mass. Events from 
the Z peak region therefore are shifted towards smaller values of m{^'^^~), thus reducing 
the cross section in and above the peak region, and increasing the rate below the Z pole. 
Due to the log(s/m|) factor, the effect of the corrections is larger in the electron case. The 
lowest order cross section is almost indistinguishable for the two lepton flavors. 

The size of the QED corrections to lepton pair production at the Tevatron becomes more 
apparent in Fig. ^ where we display the ratio of the 0{a^) and the Born cross section as a 
function of the lepton pair invariant mass. For 40 GeV < m{l^l^) < 110 GeV, the cross 
section ratio is seen to vary rapidly. Below the Z peak, QED corrections enhance the cross 
section by up to a factor 2.7 (1.9) for electrons (muons). The maximum enhancement of the 
cross section occurs at m(£+£-) ^ 75 GeV. For m{t+t) < 40 GeV (m(£+£-) > 130 GeV), 
0{a) QED corrections uniformly reduce the differential cross section by about 7% (12%) 
in the electron case, and 2.5% (~ 7%) in the muon case. Integrating over the full di- 
lepton invariant mass region, the large positive and negative corrections below and above 



Mz cancel [^]. For 40 GeV< m{i'^i ) < Mz, a large fraction of events contains a photon 

12 



with energy > 1 GeV. As we have stated before, the dominant QED radiative corrections 
are proportional to log(s/m^). The pp —>■ fi^n" cross section is therefore less affected by 
radiative corrections than the pp —>■ e+e^ rate. 

It should be emphasized that the enhanced significance of the 0{a) QED corrections 
below the Z peak is a direct consequence of the Breit-Wigner resonance of the Z boson. 
The O(a^) radiative corrections therefore should be a factor 0{{a/'n')\og{s/mf)) smaller 
than the 0{a) corrections. The effect of higher order QED corrections on the m{i~^i~) 
distribution can be estimated using the fragmentation function approach of Ref. 



this approach, the radiatively corrected cross section is obtained by convoluting the lowest 
order di-lepton cross section with a radiator function, which to all order sums the dominant 
and non-dominant logarithmic terms. Figure ^ displays the ratio of the 0{a^) cross section 
and the cross section in the fragmentation function approach as a function of m(£+£^). Only 
final state corrections are taken into account in the fragmentation function approach. As 
for the 0{a) corrections, initial state radiation contributions are expected to be small and, 
therefore, are ignored. Figure ^ shows that higher order final state QED corrections reduce 
the effect of the 0{a) corrections and are indeed of the size naively expected. In the Z peak 
region, the higher order final state corrections vary rapidly with ■m{i~^i~) and change the 
differential cross section by up to 10% (3%) in the electron (muon) case. 

In Fig. we compare the impact of the full 0{a) QED corrections (solid line) on the 
muon pair invariant mass spectrum with that of final state (dashed line) and initial state 
radiative corrections (dotted line) only. Qualitatively similar results are obtained in the 
electron case. Final state radiative corrections are seen to completely dominate over the 
entire mass range considered. They are responsible for the strong modification of the di- 
lepton invariant mass distribution. In contrast, initial state corrections are uniform and 
small (^ +0.4%). 

At small di-lepton invariant masses, photon exchange dominates and the initial - final 
state interference terms are almost completely antisymmetric in cos 6* P2|,|5B|,|55| , where 6* 



is the lepton scattering angle in the parton center of mass frame. The contribution of these 
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interference terms to the di-lepton invariant mass distribution is extremely small (0.01% - 
0.1%) for m{£~^i^) < Mz- For values of m(£+£^) sufficiently above the Z mass, initial - 
final state interference terms reduce the 0{a^) cross section by about 1%. 

Next-to- leading order QCD corrections to lepton pair production in pp collisions at Teva- 



are initial state corrections, the NLO QCD to leading order cross section ratio varies only 
slowly with the di-lepton invariant mass, similar to what we found for initial state QED 
corrections. Comparing the size of the 0{a) QED and 0{as) QCD corrections, one ob- 
serves that they are of similar magnitude above the Z peak, but have opposite sign. In 
the invariant mass range between 50 GeV and 100 GeV, QED corrections are significantly 
larger than those induced by the strong interactions. The relative importance of the QED 
corrections is due to the combined effect of mass singular logarithms associated with final 
state photon radiation, and the Z boson Breit-Wigner resonance. 

Since QED corrections strongly affect the shape of the lepton pair invariant mass distri- 
bution below the Z peak, one expects that they may also have a significant impact on other 
observables in this region. In Fig. |^, we show the forward backward asymmetry, Aps, as 
a function of the lepton pair invariant mass in the Born approximation (dashed line), and 
including 0{a) QED corrections for electron (solid line) and muon final states (dotted line). 
Here, Aps is defined by 



tron energies are known |^ to enhance the cross section by about 16% - 25%. Since these 



A 



F-B 



(13) 



FB — 



F + B 



where 




(14) 



cos 6'* is given by |lT8i|39 | 




2 



(15) 



with 
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P^ = ^{E±p.), (16) 

where E is the energy and pz is the longitudinal component of the momentum vector. In this 
definition of cos^*, the polar axis is taken to be the bisector of the proton beam momentum 
and the negative of the anti-proton beam momentum when they are boosted into the 
rest frame. In pp collisions at Tevatron energies, the fiight direction of the incoming quark 
coincides with the proton beam direction for a large fraction of the events. The definition 
of cos^* in Eq. ([T5|) has the advantage of minimizing the effects of the QCD corrections 
(see below). In the limit of vanishing di-lepton px, 0* coincides with the angle between the 
lepton and the incoming proton in the rest frame. 

Our result for ApB in the Born approximation agrees with that presented in Ref. |T6[ . 



As expected, the 0{a) QED corrections to ApB are large in the region below the Z peak. 
Since events from the Z peak, where ApB is positive and small, are shifted towards smaller 
values of m (£+£") by photon radiation, the forward backward asymmetry is significantly 
reduced in magnitude by radiative corrections for 50 GeV < m{i^i~) < 90 GeV. 

The forward backward asymmetry in the Born approximation is small at low di-lepton 
masses, because of the dominance of photon exchange and the vectorlike coupling of the 
photon to leptons. For di-lepton masses below 40 GeV, the 0{a) initial - final state in- 
terference correction terms are almost completely antisymmetric in cos^* and comprise the 
most important component of the QED corrections to ApB- In this region, the 0{a) QED 
corrections to Aps are therefore large. Initial - final state interference terms do not contain 
any mass singular contributions. As a result, the forward backward asymmetries for electron 
and muon final states are similar for m{i^i~) < 40 GeV. Details of the asymmetry in the 
low di-lepton mass region are shown in the inset of Fig. |^. Effects from purely weak cor- 
rections are not included in our calculation. They could have a non-negligible effect on the 
forward backward asymmetry at low di-lepton masses, similar to the situation encountered 
in e+e- ^ fi+fi^ @]. 

In contrast to the lepton pair invariant mass distribution, QED corrections to Aps are 
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small for m{i^i~) > 120 GeV. They reduce the forward backward asymmetry by about 1% 
in this region. Initial and final state corrections to ApB are of similar size for lepton pair 
invariant masses above the Z peak. 

Recently, the CDF Collaboration has presented a first measurement of the integrated 
forward backward asymmetry in pp e~^e~X at the Tevatron for m(e^e~) > 105 GeV, 
together with a more refined measurement in the Z peak region (75 GeV < m(e+e^) < 
105 GeV) |18|. In Table |l|, we list the experimental values, together with the theoretical 
prediction with and without 0{a) QED corrections. QED corrections are seen to increase 
the asymmetry by about 8% in the peak region. In the muon channel, the increase in ApB 
for 75 GeV < m{fi^fi^) < 105 GeV due to radiative corrections is approximately 4%. 

In the Z peak region, Aps provides a tool to measure sin^ 6'^^^ ||12|. For 75 GeV < 
m(£+£^) < 105 GeV and ^/s = 1.8 TeV, the forward backward asymmetry can to a very 
good approximation be parameterized by [|16| 



ApB = b (a - sin^ 



(17) 



both in the Born approximation and including 0{a) QED corrections. For the parameters 
a and b we find in the Born approximation 



^Born _ Q_2454, 



b^"'"" = 3.6 



(18) 



for e~^e as well as final states, and 



^O(a^) ^ ^Born ^ ^^QED ^0(^3) ^ ^Born ^ ^^QED 



(19) 



with 



Aa'^™ ^ 0.0010, 



A6^™ ^ 



(20) 



for pp e^e (7), and 



Aa^™ ^ 0.0006, 



AfeQ™ ^ -0.3 



(21) 
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for pp jj,'^ jj, (7). The change of the effective weak mixing angle due to QED radiative 
corrections is a factor 3 to 4 larger than the current experimental uncertainty, 5 sin^ 6'^^^ = 
0.00024 im. 



In Ref. |]T2|], the approximation used to estimate the electroweak corrections to ApB 
resulted in a significant dependence of the correction to sin^ O^^Jj on the infrared cutoff 
used in the calculation. In contrast, as explained in detail in Sec. II, our results are cutoff 
independent. This will make it possible to substantially reduce the theoretical uncertainty 
of the weak mixing angle extracted from future measurements of Aps at the Tevatron. 

B. Aspects of Experimental Lepton Identification and QED Radiative Corrections 



It is well-known ||3^ that the mass singular logarithmic terms which appear in higher 
orders of perturbation theory are eliminated when inclusive observables are considered. As 
explained below, the finite resolution of detectors prevents fully exclusive measurements. 
Detector effects, which we have completely ignored so far, therefore may significantly modify 
the effect of QED radiative corrections. To simulate detector acceptance, we impose the 
following transverse momentum {pj) and pseudo-rapidity [r]) cuts: 

electrons muons 
Pj.{e) > 20 GeV > 25 GeV 

\v{e)\ < 2.4 |7?(/x)| < 1.0 

In addition, we require that at least one electron (muon) is in the central part of the detector: 
\ri{e)\ < 1.1 {\ri{fi)\ < 0.6). These cuts approximately model the acceptance of the CDF 
detector for electrons and muons. Uncertainties in the energy measurements of the charged 
leptons in the detector are simulated in the calculation by Gaussian smearing of the particle 
four-momentum vector with standard deviation a which depends on the particle type and 
the detector. The numerical results presented here were calculated using a values based 
on the CDF |^ specifications. Similar results are obtained if the acceptances and energy 



resolutions of the D0 detector are used 0. 
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The granularity of the detectors and the size of the electromagnetic showers in the 
calorimeter make it difficult to discriminate between electrons and photons with a small 
opening angle. We therefore recombine the four-momentum vectors of the electron and 
photon to an effective electron four-momentum vector if both traverse the same calorimeter 
cell, assuming a calorimeter segmentation of Arj x A0 = 0.1 x 15° (0 is the azimuthal angle 
in the transverse plane). This procedure is similar to that used by the CDF Collaboration. 
The segmentation chosen corresponds to that of the central part of the CDF calorimeter . 
The D0 Collaboration uses a slightly different recombination procedure where the electron 
and photon four-momentum vectors are combined if their separation in the pseudorapidity 
- azimuthal angle plane, Ai?(e, 7) = \J {Ar]{e, 7))^ + (A0(e, 7))^, is smaller than a critical 
value, Rc- For Rc = 0.2 the numerical results obtained are similar to those found with 
the calorimeter segmentation we use (see above). 

Muons are identified in a hadron collider detector by hits in the muon chambers. In ad- 
dition, one requires that the associated track is consistent with a minimum ionizing particle. 
This limits the energy of a photon which traverses the same calorimeter cell as the muon to 
be smaller than a critical value E^. In the subsequent discussion, we assume = 2 GeV [§. 

In Fig. ^a (Fig. |^) we show how detector effects change the ratio of the 0{a^) to leading 
order differential cross sections as a function of the e~^e^ (fi^fi^) invariant mass. The finite 
energy resolution and the acceptance cuts have only a small effect on the cross section ratio. 
The lepton identification criteria, on the other hand, are found to have a large impact. 
Recombining the electron and photon four-momentum vectors if they traverse the same 
calorimeter cell greatly reduces the effect of the mass singular logarithmic terms. These 
terms survive only in the rare case when both particles are almost collinear, but hit different 
calorimeter cellsQ Although the recombination of the electron and photon momenta reduces 



^In the case where the four-momentum vectors of the two particles are recombined for AR{e, 7) < 
Rc, the mass singular terms are entirely eliminated, and the lepton mass in the logarithmic terms 
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effect of the 0{a) QED corrections, the remaining corrections are still sizeable. Below (at) 
the Z peak, they enhance (suppress) the lowest order differential cross section by up to a 
factor 1.6 (0.9) [see Fig. For m(e^e^) ^ Mz, the magnitude of the QED corrections is 
reduced from approximately 12% to 5%. 

For muon final states (see Fig. ^), the requirement of < = 2 GeV for a photon 
which traverses the same calorimeter cell as the muon reduces the hard photon part of the 
0{a^) ^^fi~{'~f) cross section. As a result, the magnitude of the QED corrections below the 
Z peak is reduced. At the Z pole the corrections remain unchanged, and for yU^/x^ masses 
larger than Mz they become more pronounced. For mlfi^fi^) > 120 GeV, QED corrections 
reduce the cross section by 12% to 14%. 

We would like to emphasize that the survival of mass singular terms in certain cases 



does not contradict the KLN theorem [^. The KLN theorem requires that mass singular 
logarithmic terms which appear in higher orders of perturbation theory are eliminated when 
inclusive observables are considered. Recombining the lepton and photon momenta for small 
opening angles an inclusive quantity is formed, and the mass singular logarithmic terms are 
eliminated in the reconstructed invariant mass distribution. On the other hand, if the 
lepton and photon momenta are not combined, one performs an exclusive measurement, 
the KLN theorem does not apply, and logarithmic terms remain present in the measured 
di-lepton invariant mass distribution. 

It should be noted that the differential cross section ratio shown in Fig. ^ becomes ill 
defined in the threshold region m (£"*"£") ^ 2p^*(£), where p^^{i) is the charged lepton 
threshold. For m (£"*"£") < 2p^*(£), the Born cross section vanishes, and the cross section 
ratio is undefined. The 0{a^) cross section is small, but non-zero, in this region. The largest 
contribution to the cross section for m{i^i~) < 2p^*(£) originates from initial state radiation 
configurations, where the leptons have a small relative opening angle and are balanced by 



is replaced by the minimum 67 invariant mass. 
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a high pt photon in the opposite hemisphere. Close to the threshold, m{i~^i~) ^ 2p™*(£), 
large logarithmic corrections are present, and for an accurate prediction of the cross section 
those corrections need to be resummed. The results of Fig. ^ in this region should therefore 
be interpreted with caution. Similar conclusions can also be drawn for the forward backward 
asymmetry in the threshold region. 

In Fig. ^ we show how detector effects affect the forward backward asymmetry for 
electron (Fig. ^) and muon final states (Fig. §3). In addition to the cuts listed at the 
beginning of this subsection, we require 



|cosr|<0.8. (22) 

For comparison, we also show the asymmetry in the Born approximation without taking any 
detector related effects into account (dotted line). The finite lepton rapidity coverage and 
the I cos 6* I cut significantly reduce the forward backward asymmetry in magnitude. Energy 
and momentum resolution effects broaden the Z peak and thus introduce a characteristic S 
type bending in ApB at m{i'^i~) ^ Mz- Analogous to the di- lepton invariant mass distri- 
bution, lepton identification requirements substantially reduce the impact of QED radiative 
corrections on the forward backward asymmetry below the Z peak. For ■m{i'^i~) > 100 GeV, 
they have only a small effect on ApB, similar to the case where no detector effects are taken 
into account. 

Although QED corrections to the forward backward asymmetry are reduced in magnitude 
for m{£^i^) < Mz by experimental lepton detection and identification requirements, they 
are still considerably larger than the NLO QCD corrections in this region. This is demon- 
strated in Fig. |l^ for the electron final state. Similar results are obtained for pp fi^^^ {'-/). 
The 0{as) QCD corrections to pp — > Z, 7* — > i^i~X are calculated in the MS scheme us- 
ing the Monte Carlo approach of Ref. |19[. The calculation generalizes that of Ref. |^ 



to include finite Z width effects and virtual photon exchange diagrams. The QCD cor- 
rections to ApB P3[ are found to be quite small. Below (above) the Z peak, the mag- 



nitude of the forward backward asymmetry is reduced by typically 6Afb/Afb ~ —0.05 
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{5Afb/Afb ~ -0.02). For 75 GeV < m(£+£-) < 105 GeV, NLO QCD corrections de- 
crease the integrated asymmetry by 5Apb/Afb ~ —0.03. QED and QCD corrections to the 
integrated forward backward asymmetry in the Z peak region have opposite signs. 

To reduce the background from heavy flavor production processes, the leptons in Z 
boson events are often required to be isolated. A lepton isolation cut typically requires the 
transverse energy in a cone of size Rq about the direction of the lepton, -E^", to be less than 
a fraction, e^;, of the lepton transverse energy -Et(^), ^-e. 

Er{i) ^ ^^^^ 

Sometimes the energy, instead of the transverse energy is used in the isolation require- 
ment, Eq. (pSl). The isolation requirement and the cut imposed on the photon energy in the 
muon case have similar effects. In Fig. |Tl|, we show how the lepton isolation requirement 
of Eq. ( p3D with Rq = 0.4 and eE = 0.1 modifies the effect of the 0{a) QED corrections 
on the di-lepton invariant mass distribution in the Z peak region. The isolation cut is seen 
to mostly affect the mass region below Mz, reducing the maximum enhancement of the 
differential cross section by QED radiative corrections from a factor ~ 1.6 to 1.2 - 1.3. In 
our calculation, for electrons, the isolation requirement is only imposed if the electron and 
photon are not recombined. 0{aas) corrections to di-lepton production are not included 
in the results presented. These corrections are expected to increase somewhat, and 
therefore will modify the effect of the isolation cut. 

In the past, the measurement of the W and Z boson cross sections has provided a test 
of perturbative QCD [^^ -^. With the large data set accumulated in the 1994-95 Tevatron 



collider run, the uncertainty associated with the integrated luminosity (~ 3.6% |^6|) becomes 
a limiting factor in this measurement. This suggests to use the measured W and Z boson 
cross sections to determine the integrated luminosity in future experiments P6| , ^ . In order 
to accurately measure the integrated luminosity, it will be necessary not only to take the 
corrections to the W and Z boson cross sections into account, but also to correct for 
higher order QED effects. 
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Experimentally, the Z boson cross section is extracted from the di-lepton cross section in 
a specified invariant mass interval around the Z boson mass, correcting for photon exchange 
and interference effects. The size of the 0[(y) QED corrections to the total di-lepton 
cross section is sensitive to the lepton identification criteria, the acceptance cuts and the 
range of the di-lepton invariant masses selected (see Fig. |TT]). In Table |I| we list the cross 
section ratio ("QED i^-factor") 

^QED ^ /24) 

0-Born ^ ' 

for 75 GeV < m(£"'"^~) < 105 GeV [i = e, /i). For comparison, we also tabulate the 
corresponding QCD K-factor, 

= ^ • (25) 

One observes that the effect of the large QED corrections found in da/dm{£'^i~) is strongly 
reduced when integrating over a range in invariant mass which is approximately centered at 
Mz- Nevertheless, the QED corrections usually are not negligible when compared with the 
0{as) QCD corrections. QCD corrections enhance the Z boson production rate, whereas 
QED effects decrease the cross section for the invariant mass window chosen here. The 
total pp e'^e'X {pp — > fi^fi^X) cross section is reduced by about 7% (3%) by QED 
radiative corrections. As we have noted before, the dominant QED correction terms are 
proportional to log(s/m^) in absence of detector related effects. Without detector effects 
taken into account, QED corrections to pp — > e^e^X thus are larger than for di-muon 
production. The recombination of electron and photon momenta when the opening angle 
between the two particles is small strongly reduces the effect of the QED corrections to 
the integrated e'^e^ cross section. In the muon case, lepton identification requirements 
increase the magnitude of the QED corrections, and they almost compensate the cross 
section enhancement originating from 0{as) QCD corrections. Requiring the lepton to be 
isolated reduces the hard photon contribution to the 0{a^) cross section, and hence increases 
the effect of the QED corrections. QCD corrections are only slightly modified by detector 
effects. 
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Since the 0{a) QED corrections and the 0{as) QCD corrections are of similar magnitude 
in the muon case when reahstic experimental conditions are taken into account, one expects 
that the 0{aas) and Ola"^) corrections are also of similar size in this channel. The 0{aas) 
corrections may thus be non- negligible in a precise determination of the integrated luminosity 
from the Z ^^fi^ cross section. 

Finite detector acceptance cuts do not significantly modify the QED corrections to 
da / dm{i~^ i~) and ApB, except in the threshold region, m{i~^i~) ^ 2p^^{i). The effect 
of the cuts can be more pronounced in other distributions. As an example, we show the ra- 
tio of the lepton transverse momentum distribution at 0{a^) and in the Born approximation 
in Fig. |T2|. All criteria which are necessary to simulate lepton detection and identification, 
except the isolation cut of Eq. (p3D, are imposed in this figure. For the CDF inspired pseu- 
dorapidity and px cuts we use in the muon case, the 2^2 phase space becomes much 
more restricted than the 2—^3 phase space close to the pt threshold. As a result, the 
cross section ratio exhibits a bump located at Pt{.Ij) ~ 30 GeV (dashed line). Replacing the 
acceptance cuts by those used for electrons, the bump in the cross section disappears (dotted 
line). For prif^) > 40 GeV, the size of the radiative corrections is almost independent of the 
pseudorapidity and transverse momentum cuts imposed. Radiative corrections smear out 
the Jacobian peak, causing a characteristic dip in the cross section ratio at pri^i) ~ Mzjl- 
However, in this region, the cross section is subject to large QCD corrections |^ which are 
not taken into account in our calculation. 

The QED corrections to the muon transverse momentum distribution reduce the cross 
section by 10 - 15% over most of the pt range. For comparison, we also display the ratio of 
differential cross sections for pp e^e~(7) in Fig. 0. Here, the 0[ol) QED corrections are 
of (9(1%), except for the Jacobian peak region, pxie) ~ 45 GeV, where they reduce the cross 
section by up to 7%. The pronounced difference in radiative corrections between electrons 
and muons is largely due to the different lepton identification requirements discussed earlier 
in this subsection. 
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C. Radiative Corrections and the Z Boson Mass 



As we have seen, final state bremsstrahlung severely distorts the Breit-Wigner shape of 
the Z resonance curve. As a result, QED corrections must be included when the Z boson 
mass is extracted from data, otherwise the mass extracted is shifted to a lower value. In 
absence of detector effects, the Z mass shift is approximately given by 



AMz ^ Tz, (26) 

o 

with /3 defined in Eq. (|T|). For Z e+e" (Z ^ AM^ ^ -110 MeV (AM^ ^ 

—60 MeV). However, as it is clear from the previous section, detector effects significantly 
modify AM^. 

The Z boson mass extracted from Tevatron experiments serves as a reference point 
when compared with the precise measurement performed at LEP. It helps to calibrate the 
electromagnetic energy scale, and to determine the electron energy resolution as well as the 
muon momentum resolution which are important for the measurement of the W mass. 

In the approximate treatment of the QED corrections used so far by the Tevatron ex- 
periments, only final state corrections are taken into account. In addition, the effects of soft 
and virtual corrections are estimated from the inclusive 0{o?^ Z —>■ i^i~{'~f) width |^ and 
the hard photon bremsstrahlung contribution W3 . 



We now study the differences in the Z boson masses extracted using the approximation 
currently employed in the experimental analysis and our complete 0{a^) QED calculation, 
and investigate the effect of the initial state radiative corrections on the Z mass shift. To 
extract the Z boson mass, we use a log-likelihood fit to the shape of the di-lepton invariant 
mass distribution in the range 81 GeV< m(£"'"£~) < 101 GeV. The templates for the m{i~^i~) 
distributions are calculated using the lowest order differential cross section, varying Mz 
between 90.6 GeV and 91.5 GeV in steps of 100 MeV. Detector effects are simulated as 
described in Sec. IIIB. No isolation cut [Eq. p3|)] is imposed on the charged leptons. The 
soft and coUinear cutoff parameters are chosen to be Sg = 10~^ and 6c = 3x 10~^. In order to 
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be able to determine AM^, it is necessary to properly include the radiation of photons with 
an energy which is of the same order as the shift in M^, using the full 2^3 phase space. 
6s and 6c, therefore, have to be smaller than about 2 x 10"^, otherwise a non-negligible 
dependence of the Z boson mass shift, AMz, on these parameters remains. 

The error on the Z mass resulting from the statistical uncertainties in the Monte Carlo 
event samples and the finite step size in varying Mz in the templates is approximately 5 MeV 
in our simulation. This is adequate for the semi-quantitative analysis reported here. It is 
straightforward to reduce the uncertainty by increasing the number of events generated and 
the number of templates used, given sufficient computing power. 

For definiteness, we concentrate on the electron channel. Results similar to those which 
we obtain are expected in the muon case. The shift in Mz induced by the QED corrections 
is determined by comparing the shape of the 0{a^) e^e~ invariant mass distribution for 
the nominal value of Mz = 91.187 GeV with that of the templates, and calculating the log- 
likelihood as a function of the Z boson mass used as input in the templates. Repeating this 
procedure 1000 times with 10,000 events each, the difference between the average of the mass 
which maximizes the log-likelihood and the nominal Z boson mass is then identified with 
the shift induced by the QED corrections. The same procedure is carried out to compute the 
Z mass shift if the approximate calculation of Ref. [l^ is used. The Z boson mass obtained 



from the complete 0{a^) cross section is found to be about 10 MeV smaller than that 
obtained using the approximate calculation. Most of the change can be attributed to the 
different treatment of the final state soft and virtual corrections in the two calculations. A 
change of 10 MeV in Mz translates into a shift of several MeV in M^ through the dependence 
of the energy scale and the momentum resolution on the Z boson mass measured [§,0. For 
the current level of precision, this small shift is unimportant. However, it cannot be ignored 
for a high-precision measurement of M^. 

In order to estimate how initial state corrections and initial - final state interference 
correction terms affect the Z boson mass, we compare the mass obtained using the full 
0{a) corrections with that extracted when final state radiative corrections are taken into 

25 



account only. The two values of Mz are found to agree within the numerical accuracy of our 
simulation. Initial state radiative corrections and initial - final state interference correction 
terms therefore contribute very little to the Z boson mass shift. As we have discussed in 
Sec. II, current fits to the PDF's do not include QED effects. This introduces theoretical 
uncertainties, such as a strong dependence of the initial state corrections on the factorization 
scheme used. However, since initial state corrections essentially do not contribute to the Z 
boson mass shift, these uncertainties will have no significant effect on the Z boson mass 
extracted. This conjecture is supported by the fact that the numerical values for the mass 
shifts in the QED MS and DIS scheme are the same. 

The Z boson mass extracted from the fit to the di-lepton invariant mass distribution also 
depends on the PDF uncertainties, and the choice of the renormalization and factorization 
scale. At present, PDF's which take into account uncertainties in their fit are not generally 
available^. We therefore only consider the scale dependence here. Changing = /i^ = 
^QED = ^QCD fro^ = s to Q"^ = 100 s decreases the fitted Z mass by 10 MeV both in 
the Born approximation and when 0{a) corrections are taken into account. This indicates 
that the Z boson mass shift caused by QED corrections is insensitive to the choice of Q^. 
The scale dependence of the fitted Z mass is eliminated when 0{as) QCD corrections are 
taken into account. 



IV. THE FORWARD BACKWARD ASYMMETRY AT THE LHC 

As we have mentioned in the Introduction, one can use sin^ ^^j^* together with rritop to 
constrain the Higgs boson mass. At LEP, sin^ ^^j^* has been measured with an accuracy 
of approximately ±0.00024 [|l|. In order to extract the Higgs boson mass with a precision 
of 6Mh/Mh ~ 30% or better, the uncertainty in svd? O^^Jj has to be reduced by at least a 
factor two. At the LHC {pp collisions at ^/s = 14 TeV), the Z i^i" cross section is 



^An approach to extract PDF's including systematic errors has recently been described in Ref. |5C 
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approximately 1.6 nb for each lepton flavor. For the projected yearly integrated luminosity 
of 100 fb~^, this results in a very large number of Z ^ events which, in principle, can 
be used to measure the forward backward asymmetry and thus sin^ 6'^^^ with extremely high 
precision [|ri|] . In this Section, we investigate the prospects to measure sin^ O^^J^ using the 



forward backward asymmetry at the LHC, taking into account the 0{a) QED and 0{as) 
QCD corrections. At LHC luminosities, it is easier to trigger on than on e+e^ pairs 

in the Z mass region ||5l| , |52[| . In our analysis, we therefore concentrate on the Z fi'^fi~ 



channel; qualitatively similar results are obtained for the electron channel. 

In pp collisions, the quark direction in the initial state has to be extracted from the boost 
direction of the di-lepton system with respect to the beam axis . The cosine of the angle 
between the lepton and the quark in the rest frame is then approximated by 

b.(£+r)| 2 



cosr 



p+{r)p-{r)-p-{r)p+{r) . (27) 



m{£+£-)^m^ {£+£-) +pU£+i-) 

For the deflnition of cos 6'* given in Eq. (^), ApB = ^ for pp collisions. 

At the LHC, the sea - sea quark flux is much larger than at the Tevatron. As a result, 
the probability, fq, that the quark direction and the boost direction of the di-lepton system 
coincide is signiflcantly smaller than one. The forward backward asymmetry is therefore 
smaller than at the Tevatron. Events with a large rapidity of the di-lepton system, y{£'^£~), 
originate from collisions where at least one of the partons carries a large fraction x of the 
proton momentum. Since valence quarks dominate at high values of x, a cut on the di-lepton 
rapidity increases fq, and thus the asymmetry and the sensitivity to the effective weak 
mixing angle. 

The forward backward asymmetry at the LHC, using Eq. (|27|) to deflne cos 6'*, and 
imposing a 

|l/(/iV)l>l (28) 
cut, is shown in Fig. O for values of m(/i+yU^) up to 250 GeV. No other cuts besides 



the y{fJ'~^fi ) cut have been imposed in Fig. |l^. Without the cut of Eq. (|28|), Aps would 
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be approximately a factor 1.25 smaller. Although the di-lepton rapidity cut enhances the 
asymmetry, it is about a factor 1.5 smaller than at the Tevatron. 

Qualitatively, the behaviour of the forward backward asymmetry as a function of the 
di-lepton invariant mass is similar to that in pp collisions. Furthermore, QED and QCD 
corrections are seen to have a quantitatively similar effect on ApB as in pp collisions. In the 
Z peak region, 75 GeV < m{^^fi^) < 105 GeV, the integrated forward backward asymmetry 
can again be parameterized by Eqs. (|T7|) and (|19|) with 

^Born _ 0.2458, 6^°"^"^ = 2.19, (29) 

AqQ™ = 0.0008, Afe'^ED _ _o.o9. (30) 

and 

^^QCD _ _o.ooil, AfeQ^^ = 0.06. (31) 



From Eqs. (|29|) and ([TsD we observe that the parameter a is essentially the same as at 
Tevatron energies, b, on the other hand, which controls the sensitivity to the weak mixing 
angle, is significantly reduced. QED corrections increase the integrated asymmetry in the 
peak region by about 6Afb/Afb ~ 0.02, and slightly reduce the sensitivity to sin^6'gj^*. 
QCD corrections are found to reduce it by approximately 6Afb/Afb ~ —0.05. QED and 
QCD corrections to the integrated forward backward asymmetry in the Z peak region have 
opposite sign, as at the Tevatron. 

Using Eqs. (|29| ) and (|30| ) together with Eqs. ( |T^ ) and (0), it is now straightforward 
to estimate the error expected for sin^ 6'^^^ from a measurement of the forward backward 
asymmetry in the Z peak region at the LHC. For an integrated luminosity of 100 fb~^, we 
find that it should be possible to measure sin^ ^'j^* with a statistical precision of 

6 sin^ 9[ff = 3.9 xlO-^ (32) 

Both, NLO QCD and QED corrections have been taken into account in this estimate. The 
\y{fi~^fi~')\ > 1 cut improves the precision for sin^ 6*^^^* by about 10%. Our result is about 35% 
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better than the estimate given in Ref. |TT|. The shift in ApB introduced by the combined 
QED and QCD radiative corrections is about a factor 7 larger than the expected statistical 
error. 

The estimate of Eq. ( p^ ) has been obtained assuming full rapidity coverage for the 
muons. The proposed FELIX experiment |Q is expected to achieve this. However, FELIX 



will operate at a reduced luminosity of at most C = 10^^ cm~^s^^, corresponding to a yearly 
integrated luminosity of 10 fb~^ at best. For 10 fb~^ the expected precision is 5sm^9^^J^^ ~ 
1.2 X 10~^. In both the ATLAS and CMS detector, muons can only be detected for pseudo- 



rapidities < 2.4 pll , |52| . In Fig. |Tj we display the forward backward asymmetry at the 

LHC imposing a \ri{^)\ < 2.4 cut in addition to the di-lepton rapidity cut of Eq. (pSf). The 
finite rapidity range covered by the detector is seen to dramatically reduce the asymmetry. 
In the region around the Z pole, the integrated forward backward asymmetry is again an 
approximately linear function of sin^ (see Eq. (p!7|)) with 

^Born _ 0.2464, 6^°™ = 0.72. (33) 

QED and QCD radiative corrections shift these values by 

Aa^™ = 0.0024, Afo^^^ = -0.07, (34) 

and 

Aa<^^° = 0.0067, Aft^™ = -0.27, (35) 

respectively. The parameter 6, which directly controls the sensitivity to sin^ ^e//5 is reduced 
by about a factor 3 by the finite rapidity acceptance. QED corrections further reduce h 
by approximately 10%, and QCD corrections by an additional 30%. The finite rapidity 
coverage also results in a reduction of the total Z boson cross section by roughly a factor 5. 
As a result, the uncertainty expected for sin^ O^^J^^ with 100 fb~^ increases by more than a 
factor 10 to 

(5sin2 e^^ff = 4.4 X 10"^ for \r]{fi)\ < 2.4. (36) 
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The shift in ApB introduced by the combined QCD and QED radiative corrections is about 
a factor 1.5 larger than the statistical error expected. 

The rapidity range covered by the electromagnetic calorimeter and the tracking system 



is very similar to that of the muon system For e^e~ production, one therefore does 

not expect to measure sin^ ^^j^* with a higher precision than in the muon channel. As in the 
Tevatron case discussed in Sec. IIIA, the QED corrections to Afb in absence of detector 
effects are more pronounced in the electron case. A pt(^) > 20 GeV cut has essentially no 
effect on the forward backward asymmetry in the Z peak region. 

The precision expected for sin^ 6'^^^* from LHC experiments should be compared with the 
accuracy from current LEP and SLC data and with the sensitivity expected from future 
experiments at the SLC and the Tevatron. The combined uncertainty of sin^ O^^jj- from LEP 



and SLC experiments is approximately 2.4 x 10 . With the planned luminosity upgrade [55 



one hopes to collect 3 x 10^ Z boson events at the SLC. This would allow to measure sin^ ^^^^ 
from the left right asymmetry with a precision of about 1.2 x 10~^, which is similar to the one 
attainable by FELIX with 10 fb^^. At the Tevatron, with the same integrated luminosity, 
one expects an uncertainty of 2.3 x 10~^ for sin^ 6''^j^ ^ per experiment. In order to improve 
the precision beyond that expected from future SLC and Tevatron experiments, it will be 
necessary to detect leptons in the very forward pseudorapidity range, \ri\ = 3.0 — 5.0 at the 
LHC when it operates at the design luminosity of C = 10^^ cm~^s~^. 



V. SUMMARY AND CONCLUSIONS 



In a precision measurement of in hadronic collisions, a simultaneous determination 
of Mz in di-lepton production is required for calibration purposes. The forward backward 
asymmetry makes it possible to determine sin^ 6'^^'^ with high precision. Both measurements 
help to constrain the Higgs boson mass from radiative corrections. In order to perform these 
high precision measurements, it is crucial to fully control higher order QCD and electroweak 
corrections. In this paper we have presented a calculation of di-lepton production in hadronic 
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collisions based on a combination of analytic and Monte Carlo integration techniques which 
includes initial and final state 0{a) QED corrections. Previous calculations |T3| . |I^ have 
been based on the final state photonic corrections, estimating the virtual corrections indi- 
rectly from the inclusive O(a^) Z ^+£^(7) width and the hard photon bremsstrahlung 
contribution. 

Due to mass singular logarithmic terms associated with final state photon radiation in 
the limit where the photon is coUinear with one of the leptons, final state radiation effects 
dominate. Initial state corrections were found to be small after factorizing the corresponding 
collinear singularities into the parton distribution functions. QED corrections to the evolu- 
tion of the parton distribution functions and purely weak corrections are not included in our 
calculation; they are expected to be small. Initial state QED corrections are uniform over 
the entire di-lepton invariant mass range. In contrast, final state corrections vary rapidly 
with m{i~^i~), and strongly modify the shape of the invariant mass distribution as a large 
fraction of the events shifts from the Z boson peak to lower invariant masses (see Figs. |^ 
and Below Mz, radiative corrections enhance the cross section by up to a factor 2.7 (1.9) 
for electrons (muons). 

QED corrections also strongly reduce the magnitude of the forward backward asymmetry, 
ApB, for di-lepton invariant masses between 50 GeV and 90 GeV. In the Z peak region, 
75 GeV < m{i~^£~) < 105 GeV, they enhance the integrated forward backward asymmetry 
by up to 8%. 

When detector effects are taken into account, the effect of the mass singular logarithmic 
terms in the electron case is strongly reduced. The granularity of the detector and the size 
of the electromagnetic showers in the calorimeter make it difficult to discriminate between 
electrons and photons with a small opening angle. One therefore combines the electron and 
photon four momentum vectors if both particles traverse the same calorimeter cell. In the 
muon case, the energy of the photon is required to be smaller than a critical value, E^, if 
both particles traverse the same calorimeter cell, and mass singular terms survive. Removing 
energetic photons reduces (enhances) the effect of the 0{a) corrections below (above) Mz- 
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Detector effects are also found to considerably decrease the size of the QED corrections to 
the forward backward asymmetry below the Z peak for — > fi^fi~{'-f). 

QED corrections have a significant impact on the di-lepton cross section in the Z peak 
region, and the Z mass extracted from experiment. In future Tevatron runs, the total 
W/Z cross section may be used as a luminosity monitor [^. As shown in Table y, QED 
corrections can reduce the di-lepton cross section in the Z peak region by up to 10%. Final 
state radiative corrections are known P,^] to substantially shift the Z boson mass. The 
Z boson mass extracted from our 0{a^) invariant mass distribution was found to be 



about 10 MeV smaller than that obtained using the approximate calculation of Ref. |]13 
Initial state corrections and initial - final state interference terms only marginally influence 
the amount the Z boson mass is shifted. The contribution of the QED corrections to the 
PDF's is expected to be of the size of the initial state radiative corrections that are included 
in our calculation. It is unlikely to be a limiting factor in the determination of the Z (and 
W) boson mass in hadronic collisions. 

For the current level of precision, the approximate calculation of Ref. ||TB[ appears to be 
adequate. The small difference in the Z boson mass obtained in the complete 0{a^) and the 
approximate calculation, however, cannot be ignored if one attempts to measure the W mass 
with high precision at hadron colliders. This also raises the question of how strongly multiple 
final state photon radiation influences the measured Z boson mass. Using the fragmentation 
function approach, we have shown that higher order QED corrections non-trivially modify 
the shape of the di-lepton invariant mass distribution. They may introduce an additional 
shift of Mz by (9(10 MeV), and may have a non-negligible impact on the forward backward 



asymmetry. So far, only partial calculations exist [£6|. A more complete understanding of 
multiple photon radiation is warranted. 

Finally, we studied the forward backward asymmetry at the LHC The very large number 
of Z bosons produced at the LHC offers an opportunity to accurately measure sin^ from 
ApB- For the forward backward asymmetry to be non-zero in pp collisions, the scattering 
angle has to be defined with respect to the boost direction of the lepton pair along the 
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beam axis. Imposing a \y{£'^£~)\ > 1 cut reduces the fraction of events where the quark 
direction is misidentified. It enhances the asymmetry by a factor 1.25, and thus improves 
the sensitivity to sin e'-^Jf by about 10%. With a detector possessing full rapidity coverage 
for leptons, sin^ ^^J^ can in principle be measured with a precision of 5 sin^ elff = 3.9 X 10-5 
if an integrated luminosity of 100 fb"^ is achieved. The shift in ^4^^ introduced by QED and 
QCD radiative corrections is about one order of magnitude larger than the statistical error 
expected. The finite lepton rapidity coverage of the ATLAS and CMS detectors strongly 
reduces ApB and the number of Z bosons produced, which results in an increase of the 
uncertainty in sin^ by about a factor 10. In order to significantly improve the precision 
for sin^ ^^.J^ beyond that expected from future SLC and Tevatron experiments, it will thus 
be necessary to detect electrons and muons in the very forward pseudorapidity range, |?7| = 
3.0 — 5.0, at the LHC, and to achieve an integrated luminosity of C(100 fb"^). 
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TABLES 

TABLE L The integrated forward backward asymmetry, Aps, inpp — > e^e~X at ^/s = 1.8 TeV 
for 75 GeV < m(e^e~) < 105 GeV and m(e^e^) > 105 GeV. Shown are the SM predictions with 
and without 0{a) QED corrections together with the experimental values of Ref. [^]. The uncer- 
tainties hsted for the theoretical results represent the statistical error of the Monte Carlo integra- 
tion. 



75 GeV < m{e+e~) < 105 GeV m(e+e-) > 105 GeV 
g'" 0.048 ± 0.001 0.523 ± 0.001 

A^^^^^ 0.052 ± 0.001 0.528 ± 0.001 

A7{^- 0.070 ±0.016 0.43 ±0.10 



TABLE II. The cross section ratios K'^^^ = ct^("')/cj^°™ and K'^^^ = a^("^Vf^^°™ for 
pp £+rX {i = e, fi) at ^ = 1.8 TeV with 75 GeV < m{£+i-) < 105 GeV. Shown are the 
predictions for three cases: without taking any detector effects into account ("no detector effects"), 
with the detector effects described in the text and no lepton isolation cut ("with detector effects, 
no lepton isolation"), and finally adding lepton isolation [see Eq. (p3|)] ("with detector effects, with 
lepton isolation"). 







no detector effects 


with detector effects 
no lepton isolation with lepton isolation 


KQED _ 


e+e-X) 


0.93 


0.98 
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j^QED _ 


. ^i+A^-X) 
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KQCD (^pp _ 


£+£-X) 
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FIG. 1. The pp £+£-{-y), {£ = e, /x) cross section for = 1.8 TeV and 75 GeV < m{i+£-) 
< 105 GeV as a function of a) 6c for 5s = 0.01, and b) Sg for Sc = 0.0005, including initial state 
radiation corrections only. Shown are a(2 2) — (T(Born), a{2 3), and C7(NL0) — c7(Born). 
(t(NLO) denotes the 0{a^) cross section. 
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FIG. 2. The cross section a) a{pp — e'^e~{'y)) and b) a{pp — /j,'^ iJ,~ {'j)) as a function of (5s, in- 
cluding final state radiation corrections only, for ^/s = 1.8 TeV and 75 GeV < ■m{i~^i~) < 105 GeV. 
Shown are the 2 — > 2 and 2 — > 3 contributions, and the total 0{a^) cross section. The solid line 
represents the Born cross section. 
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FIG. 3. The lepton pair invariant mass distribution for pp — at ^/s = 1.8 TeV in the 

vicinity of the Z peak. The sohd (dotted) line shows d(j/dm{i'^i~) for electron (muon) final states 
including 0{a) QED corrections. The dashed lines gives the £'^i~ Born cross section. 



42 



3.0 




250 



m(l^l ) (GeV) 

FIG. 4. Ratio of the 0{a^) and lowest order differential cross sections as a function of the 
di-lcpton invariant mass ioi pp i^i~{'y) at ^/s = 1.8 TcV. The solid line shows the result obtained 
for final state electrons, whereas the dashed line displays the cross section ratio for muons. 
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FIG. 5. Ratio of the O(a^) cross section and the cross section obtained in the fragmenta- 
tion function approacli {(t^^) as a function of the di-lepton invariant mass for pp — > i^i^X at 
■s/s = 1.8 TeV. The solid line shows the result obtained for final state electrons, whereas the dashed 
line displays the cross section ratio for muons. In the fragmentation function approach, only final 
state corrections are taken into account. 
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FIG. 6. Ratio of the 0{a^) and lowest order differential cross sections as a function of the 
di-muon invariant mass for pp fi'^fi~{^) at ^/s = 1.8 TeV. The solid line gives the result for the 
full set of 0{a^) QED diagrams. The dashed and dotted lines show the ratio obtained taking only 
final state and initial state corrections, respectively, into account. 
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FIG. 7. The forward backward asymmetry, ApB, as a function of the di-Iepton invariant mass 
for pp — ^+€"(7) at i/i =1.8 TeV. The solid and dotted lines show the forward backward asym- 
metry including 0{a) QED corrections for electrons and muons, respectively. The dashed line 
displays the lowest order prediction of ^i^B. The inset provides a closeup of ApB the low mass 
region. 
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FIG. 8. Ratio of the 0{a^) and lowest order differential cross sections as a function of the 
di-lepton invariant mass for a) pp — > e^e~{^) and b) pp — /j,'^ fj,~ {'j) at ^/s = 1.8 TeV. The solid 
(dashed) lines show the cross section ratio with (without) the detector effects described in the text. 
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FIG. 9. The forward backward asymmetry, ApB, a) ioi pp — > e+e~(7) and b) ioi pp /U+/U~(7) 
at ^/s = 1.8 TeV as a function of the di-lepton invariant mass. The sohd hnes show the result of 
the 0{a^) calculation including detector effects (see text for details). The dashed and dotted lines 
represent the forward backward asymmetry in the Born approximation with and without detector 
effects. 
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FIG. 10. The forward backward asymmetry, ApB-, including detector effects (see text for de- 
tails) as a function of tlic e+e^ invariant mass for er X at \fs = 1.8 TeV. The curves 
are for the forward backward asymmetry in the Born approximation (dotted line), including 0(a) 
QED corrections (solid line), and including 0{a^ QCD corrections (dashed line). 
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FIG. 11. Ratio of the 0{a^) and lowest order differential cross sections, including detector 
effects (see text for details), as a function of the di-lepton invariant mass for pp — > £~^£~{'y) at 
^/s = 1.8 TeV in the Z peak region. The solid and dotted lines show the cross section ratio 
without imposing a lepton isolation cut for electrons and muons, respectively. The short-dashed 
and long-dashed lines give the result imposing in addition the isolation requirement of Eq. (23) 
with Rq = 0.4 and ce = 0.1. 
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FIG. 12. Ratio of the 0{a^) and lowest order differential cross sections, including detector 

effects (sec text for details), as a function of the Icpton transverse momentum in the reaction 
pp — i~^£~{'y) at -s/s = 1.8 TeV. The solid and dashed lines show the cross section ratio for elec- 
trons and muons, respectively, employing the acceptance cuts listed in the text. The dotted line 
displays the results for muons if the same pseudorapidity and px cuts as for electrons are used 
imim^ = max(|7?(^+)|, \r,{n\y, |7?(^)Uin = min(|r?(£+)|, \r,{e-)\)]. 
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FIG. 13. The forward backward asymmetry, ApB-, as a function of the ii^ pT invariant mass for 
/i+/i^(7) at -y/i = 14 TcV. The solid and dotted lines show the forward backward asymmetry 
including Oip) QED and 0{a.s) QCD corrections, respectively. The dashed line displays the lowest 
order prediction of ApB- A |y(|Lt"^jLt~)| > 1 cut is imposed on the rapidity of the muon pair. No 
detector effects are included here. 
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FIG. 14. The forward backward asymmetry, ApB, as a function of the invariant mass 

for pp fi'^fi^l'j) at ^/s = 14 TeV. A < 2.4 cut is imposed in addition to the |y(/7,+yU^)| > 1 

cut. The solid and dotted lines show the forward backward asymmetry including 0{a) QED and 
0{as) QCD corrections, respectively. The dashed line displays the lowest order prediction of ApB- 
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